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Ground state energies and superfluid gaps are calculated 
for degenerate Fermi systems interacting via long attractive 
scattering lengths such as cold atomic gases, neutron and nu- 
clear matter. In the intermediate region of densities, where 
the interparticle spacing (~ is longer than the range 

of the interaction but shorter than the scattering length, the 
superfluid gaps and the energy per particle are found to be 
proportional to the Fermi energy and thus differs from the 
dilute and high density limits. The attractive potential in- 
crease linearly with the spin-isospin or hyperspin statistical 
factor such that, e.g., symmetric nuclear matter undergoes 
spinodal decomposition and collapses whereas neutron mat- 
ter and Fermionic atomic gases with two hyperspin states are 
mechanically stable in the intermediate density region. The 
regions of spinodal instabilities in the resulting phase diagram 
are reduced and do not prevent a superfluid transition. 



I. INTRODUCTION 

Dilute degenerate Fermi systems with long scatter- 
ing lengths are of interest for nuclear and neutron star 
matter (see, e.g., Recently, also dilute systems 

of cold Fermionic atoms have been trapped |||. The 
number density is sufficient for degeneracy to be ob- 
served and superfluidity is expected at critical temper- 
atures similar to the onset of Bose-Einstein condensates, 
~ 10 — lOOnK. S-wave scattering lengths can be very 
long, e.g.[] a = —2160 Bohr radii for triplet 6 Li and 
a ~ —18.8 fm for neutron-neutron interactions. These 
scattering lengths \a\ are much longer than typical range 
of the potentials, respectively R ~ 1 fm for strong inter- 
actions and R ~ 10 — 100 A for van der Waals forces. 
Generally, when \a\ S> R three density regimes natu- 
rally emerges: the low density (or dilute, k^^a]), the 
high density (hp £R), and the intermediate density re- 
gion (R^kp 1 ^]). 

The latter "novel" region of densities is the object of 
study here. It will be shown that previous conjectures 
based on extrapolations from the dilute limit fail. In- 
stead it is found that both the energy per particle and 
superfluid gaps scale with the Fermi energy. They de- 
pend only on statistical factors but not on the scattering 



length, range or other details of the interaction. Con- 
sequently, phase diagrams are dramatically altered and 
the stability criteria differ so that two spin systems are 
actually stable. 

The manuscript is organized as follows. In Sec. II the 
scaling and poles of the effective scattering amplitude 
are studied in a homogeneous many-body system going 
from dilute to intermediate densities. In Sec. Ill the 
ground state energies are calculated at intermediate den- 
sities and compared to well known results from the dilute 
and high density limits. In Sec. IV extensions to finite 
temperatures are discussed and a phase diagram is con- 
structed displaying regions of superfluidity and spinodal 
instabilities. In Sec. V the properties of finite systems 
of Fermions are investigated as they are relevant for cur- 
rent experiments with magnetically trapped cold atoms. 
Finally, a summary and conclusion is given. 



II. THE EFFECTIVE SCATTERING 
AMPLITUDE AND SUPERFLUIDITY 

Consider a homogeneous many-body system of 
fermions of mass m at density: p = iskp/6ir 2 , where 
v is the statistical factor, e.g. v = 4 for symmetric nu- 
clear matter and v — 2 for neutron matter as well as 
for the 40 K atomic gas with two hyperspin states cur- 
rently studied at JILA [pj. The scattering lengths and 
Fermi momentum fcj? are assumed the same for all spin- 
isospin/hyperspin components in the system but interest- 
ing effects of varying the relative densities of the various 
components will be discussed at the end. Particles are as- 
sumed non-relativistic and to interact through attractive 
two-body contact interactions. The details of the poten- 
tial is not important, only its range ~ R and scattering 
length a. We shall be particular interested in cases where 
R <C |a|, which occur when, for example, the two-body 
potential can almost support a bound state or resonance. 

At dilute or intermediate densities the particles inter- 
act via short range interactions that appear singular on 
length scales of order the interparticle distance ~ kp 1 . 
Such systems are best described by resumming the mul- 
tiple interactions in terms of the scattering amplitude. 
The Galitskii's integral equations [pj for the effective 
two-particle interaction or scattering amplitude in the 
medium is given by the ladder resummation 

r( P , P ', P) = r ( P , P ', P) + m ]T r ( P , k, p) 



The sign convention of negative scattering length for at- 
tractive potentials is used. Also % = c = 1. 
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Here, r = Ana/m is the s-wave scattering amplitude 
in vacuum; the total energy of the pair in the center 
of mass is K 2 /m; p,k, p' are the relative momentum of 
the pair of interacting particles in the initial, interme- 
diate and final states respectively, and P the total mo- 
mentum; A(P,k) = +1 for particle-particle propagation 
(|P ± k| > kp), N(P, k) = —1 for hole-hole propagation 
(|P±k| < kp), and zero otherwise. For spin independent 
interactions the amplitudes contain a factor (1 — 5„ lV2 ) 
that takes exchange into account between identical spins 
vi = v 2 . 

For obtaining BCS gaps it is sufficient to study pairs 
with P = where Eq. (fil) is simply 



r = r + 2m 
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Note the factor of 2 due to particle-particle and hole- 
hole propagation each contributing by the same amount 
in non-dense systems. The ladder resummation implicit 
in Eq. ([!]) insures that only momenta smaller than Fermi 
momenta enter. Tq vary on momentum scales ~ 1/R 3> 
kp only and can therefore be considered constant at low 
and intermediate densities. Eq. (|J) is then easily solved 
for momenta near the Fermi surface 
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The in-medium scattering amplitude has a pole due to 
Cooper pairing when 
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where Ep — k 2 F /2m is the Fermi energy. The critical 



temperature is T c = ('j/n) A, where 7 



and C = 



0.577 is Euler's constant. Eq. (|4|) is the BCS gap in the 
dilute limit which agrees with gaps calculated in H . 

However, Gorkov et al. || showed that spin fluctua- 
tions lead to a higher order correction ~ (kpa) 2 in the 
denominator of Eq. (||) that is amplified by logarithmic 
terms ~ ln(A) ~ 1/kpa. It contributes by a (negative) 
constant in the exponent and leads to a reduction of the 
gap in the dilute limit by a factor (4e) 1//3 = 2.215... as 
compared to Eq. (Q) for two spins. Generally for v spins, 
isospins or hyperspins the gap is ||] 
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In the intermediate density region pairing must still 
occur since interaction are attractive. The validity of the 
expressions of Eqs. (f|,||) in this density regime will be 
discussed further below. They predict that in the limits 
a — > — 00 and R — > the gap cannot depend on either 
a, R or other details of the potential. For dimensional 
reasons the gap can therefore only be proportional to the 
Fermi energy. 



III. GROUND STATE ENERGIES 

The ground state energy is another crucial property of 
the system. In terms of the on-shell effective scattering 
amplitude it is 0,0 



E = H^ + \ E r(p, P ,p)(i-^ 2 ) 

^£r (P)P ,P). (6) 
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Here, N = V p is the number of particles and the summa- 
tions v\,V2 include spin and isospin or hyperspin states. 
The factor (1 - 

S Ul v 2 ) m the amplitude due to exchange 
has now been written explicitly. Antisymmetrization 
of the wave-function prevents identical particles to be 
in relative s-states. At low and intermediate densities, 
kpR <C 1, the exchange term is 1/v of the direct one for 
spin independent interactions. 

Before investigating the novel intermediate density re- 
gion, a brief review of results at low and high densities is 
given. 



A. Low densities (dilute): kp\a\ <C 1 

At low densities, kp\a\ <C 1, gaps are small and have 
little effect on the total energy of the system. Expanding 
the effective scattering amplitude of Eq. (|l|) in the small 
quantity kpa, the energy per particle is obtained from 
Eq. (^J) by summing over momenta of the two interacting 
particles 
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It consists of respectively the Fermi kinetic energy, the 
standard dilute pseudo-potential j^] proportional to the 
scattering length and density, and orders (k F a) 2 Q and 
higher @. 

The zero temperature hydrodynamic sound speed 
squared can at low temperatures be expressed as 



1 (dP\ 1 d ( 2 dE/N\ 



m \dp J m 
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(8) 



With the energy per particle of Eq. (Q) at low densities, 
the sound speed can be expanded as 



1 + -{v - l)k F a + 



(9) 



where vp = kp/m is the Fermi velocity. It is com- 
monly conjectured from the first two orders that the 
Fermion (and Bose) gases undergo spinodal instability 
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when the sound speed squared becomes negative, which 
occurs when k F a^ — ir /2{v — 1). However, at the same 
densities the dilute approximation leading to Eq. (Q) fails 
and so does the conjecture as will be shown below. 



B. High densities, kpR ~3> 1 

At high densities, k F R}tl, the particle potentials over- 
lap and each particle experience on average the volume 
integral of the potentials. The energy per nucleon con- 
sists of the Fermi kinetic energy and the Hartree-Fock 
potential, (see, e.g. @ Eq. 40.14): 
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The latter term is the exchange energy which vanishes 
at very high densities, k F R ~ R/tq 3> 1, leaving the 
Hartree potential term only. At lower densities, k F R ~ 
R/ r o <^ 1) it is identical to the first integral, i.e., the 
Hartree direct term is v times the Fock exchange term as 
also found in the dilute limit, Eqs. (0-0). 

As shown in pH , the Hartree potential is consider- 
ably less attractive than the dilute potential. In fact it 
vanishes when the range of the interaction goes to zero 
and the scattering length to infinity. Take for example a 
square well potential of range R and depth — Vq. Long 
scattering lengths requires VqR 2 — ► 7t 2 /4to, and there- 
fore the Hartree potential is oc Voi? 3 ~ R — > 0. Only in 
the Born approximation do the Hartree ( |l0| ) and dilute 
potentials (0) coincide since the Born scattering length 
is a Born = (m/47r) / d 3 rV(r). 

Short range repulsion complicate the high density 
limit. In nuclear and atomic system the repulsive core 
is only of slightly shorter range than the attractive force. 
It makes the liquid strongly correlated and the Hartree- 
Fock approximation fails [001 ■ How the short range 
repulsion turn the attraction to repulsion at these even 
higher densities will, however, not affect the intermediate 
density region. 



C. Intermediate densities, \a\ kp 1 3> R 

At intermediate densities, \a\ ^> kp 1 ^> R, the scatter- 
ing length expansion in Eq. (7|) breaks down. Brueckner, 
Bethe and Goldstone (ll]] pioneered such studies for nu- 
clear matter and 3 He where the range of interactions, 
scattering lengths and repulsive cores all are compara- 
ble in magnitude. In our case the range of interaction is 
small, kpR -C 1, and therefore all particle-hole diagrams 
are negligible. Higher order particle-particle and hole- 
hole diagrams do contribute by orders of ~ TQ(mk F T) n , 
It is evident from Eq. (^) that at intermediate densi- 
ties r no longer is proportional to Tq or the scattering 



length but instead T oc (mkp) -1 . Due to the very re- 
stricted phase space such higher order terms are usu- 
ally neglected as in standard Brueckner theory. T of 
Eq. (Q) can therefore be considered as a resummation 
of an important class of diagrams. The Cooper instabil- 
ity complicates the calculation of T. If the gap is small 
the instability occurs only for pairs near the Fermi sys- 
tem with opposite momenta and spin and the effect on 
the total energy is small. The momentum dependence of 
the effective scattering amplitude also complicates a self- 
consistent calculation. These complications can be dealt 
with by approximating T by its momentum average value 
in Eq. ([!]). The momentum integrals are then analogous 
to those in the dilute limit (m) , and one obtains from Eq. 
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This expression is valid for dilute systems, where it re- 
produces Eq. ([?]), and approximately valid within the 
Galitskii ladder resummation at intermediate densities, 
R <C kp 1 <C \a\, where it reduces to 
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(12) 



with c x = 35/9(11 - 2 In 2) ~ 0.40 and v c = 1 + 3/5ci ~ 
2.5. Both the attractive and the kinetic part of the energy 
per particle are proportional to the Fermi energy at these 
intermediate energies as found for the ga ps a bove. 

The other remarkable feature of Eq. ([12]) is that the 
energy per particle changes sign for a critical number 
of degrees of freedom, v c ~ 2.5. Fermi systems with 
more degrees of freedom such as symmetric nuclear mat- 
ter have negative energy per particle and will therefore 
implode, undergo spinodal decomposition and fragment 
fll2| . Contrarily, systems with v^v c such as neutron mat- 
ter and atomic systems with only two hyperspins have 
positive energy per particle and will therefore explode, if 
not contained. This is also evident from the sound speed 
squared which from Eqs. (|§|) and (|lj) becomes 



(5/9)ci(j/ c - v)v 



(13) 



Calculations for pure neutron matter and symmetric 
nuclear matter at low densities by variational Monte 
Carlo |lj] and in neutron matter by Pade' approximants 
to R-matrix calculations | fl3[ independently confirm the 
above results approximately in a limited range of inter- 
mediate densities. In the density range po~P^M~ 3 — 
10 _4 po the energy per particle is positive for neutron 
matter and negative for symmetric nuclear matter [0. 
They scale approximately with p 2 / 3 with coefficients 
compatible with Eq. (|l2|). In symmetric nuclear matter 
the intermediate density regime is, however, limited since 
protons also interact through the triplet channel, which 
has a shorter repulsive scattering length a t ~ 5.4 fm, be- 
sides the singlet one, a s = —18.8 fm, relevant for neutron 
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matter. Never-the-less, the ladder resummation in the 
Galitskii integral equation Eqs. ( [Ti"l]l^ ) are supported by 
dimensional arguments and quantitatively it successfully 
predicts v c between that of neutron and symmetric nu- 
clear matter. The ladder resummation therefore seems to 
include the most important class of diagrams. However, 
even small corrections can be important for the magni- 
tude of the gap because they appear in the exponent as, 
e.g., found for induced interaction (compare Eq. (g) with 
(U)). In addition, the superfluidity decrease the energy 
of the system by ~ A 2 /2Ep, which can be significant at 
intermediate densities if gaps really are as large as the 
Fermi energy. 

The interesting feature of the intermediate density re- 
gion, that energies and gaps are independent of the scat- 
tering length, leads to the remarkable fact that the sys- 
tem is also insensitive to whether the scattering length 
goes to plus or minus infinity. In other words, a many 
particle system is insensitive to whether the two-body 
system has a marginally bound state just above or below 
threshold; the two-body bound state or resonance is dis- 
solved in matter at sufficiently high density, fc^|a|^,l. For 
positive scattering lengths a pair condensate of molecules 
may form at low densities but they dissolve at intermedi- 
ate densities when the range of the two-body wavefunc- 
tion exceeds the interparticle distance. 
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FIG. 1. Ground state energy and superfluid gaps for a de- 
generate system of fermions interacting via attractive forces 
with R < |o|. The energy per particle E/N (Eq. @, full 
curves and the BCS gap A (Eq. (g), dashed curves) are 
plotted in units of the Fermi energy as function of density 
for v = 2 and v = 4. Dotted curves to the right show quali- 
tatively the gap and energy per particle at high density (see 
text) as a function of IzfR (upper axis). 

In Fig. 1 the energy per particle is shown as function 
of density extending from the dilute and intermediate 
densities, Eqs. (0|j_j|), to high densities, Eq. @. At 
low densities the Fermi kinetic energy dominates but at 



intermediate densities, R^jkp ~|a|, the attractive poten- 
tial lower the energy by an amount proportional to the 
statistical factor. The two cases v — 2 and v = 4 are 
seen to saturate at positive and negative energies respec- 
tively. In the high density limit the attractive (Hartree) 
potential of Eq. ( |Io| ) dominates and will lead to collapse 
of all Fermi systems in the absence of repulsive cores. 
Fig. 1 also shows the superfluid gaps of Eqs. 
for dilute and intermediate density Fermi systems. 
At low densities they decrease exponentially as A ~ 
Ep exp(— 2/irkp\a\) whereas at intermediate densities the 
gaps are a finite fraction of the Fermi energy. At high 
densities the gap generally decreases rapidly with density 
]l5|.[f[. For example, for an attractive square well poten- 
tial of range R and depth Vq with long scattering length 
(or marginally bound state), i.e. VqR 2 ~ 7r 2 /4m, the gap 
decrease exponentially as A ~ exp(— AkpR/ir). When 
\a\ ^ R plateaus appear for (E/N)/Ep and A/Ep. 
Since Ep also decrease with decreasing density, the gap 
itself is narrowly peaked near kp ~ l/R as found in nu- 
clear and neutron matter p5[ . 

Information on the density dependence can be ob- 
tained independently from calculations within the 
Wigner-Seitz cell approximation that has recently been 
employed for the strongly correlated nuclear liquid Jl6| . 
The periodic boundary condition is a computational con- 
venience which contains the important scale for nucleon- 
nucleon correlations given by the interparticle spacing. It 
naturally gives the correct low density Eq. (0) and high 
density Eq. (10) limits. At intermediate densities one 
obtains 
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(14) 



Finite crystal momenta complicates the calculation of Ci- 
A lower (but reasonable) estimate C2 — 0.25 can be cal- 
culated. The potential energy in Eq. (|lj) is also pro- 
portional to the kinetic one as found in Eq. ([l2]) and 
of similar magnitude. The scaling with i/ -1 / 3 arise be- 
cause energies scale with the square of the inverse particle 
spacing, r^ 2 , in the Wigner-Seitz cell approximation, and 
p = vkp/6n 2 = (47rrQ/3) _1 . The energy per particle can 
be calculated at all densities and finite values of a and 
R and the cross over between the three density regimes 
generally confirm the energy per particle shown in Fig. 
1. 



IV. PHASE DIAGRAM 

Constructing a phase diagram from the low tempera- 
ture degenerate regime to the high temperature classi- 
cal regime requires a finite temperature generalization. 
For illustration we shall follow the procedure as in Ref. 
pl[ and employ the high temperature approximation for 
the additional thermal pressure. At high temperatures 
quantal effects are negligible and the energy per particle 
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is simply given by the classical value E/N ~ 3T/2. The 
isothermal sound speed is within this approximation 



m 



(15) 



where the zero temperature sound speed is given by Eq. 
(0) with energy per particle from Eqs. (|TT|,f7|). The spin- 
odal instability condition, st = 0, determines the line of 
collapse T(p) for long wavelength density fluctuations. 
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FIG. 2. Phase diagram at low and intermediate densities 
(p = vkp /6tt) for a gas of fermions interacting via a long (at- 
tractive) scattering length a. Spinodal lines are shown with 
full curves and the superfluid transition by dashed curves for 
various number of spin states v as labeled. The area con- 
strained to the lower right corner are the spinodally unstable 
and superfluid regions. As systems with two spin states only 
are stable at intermediate densities the v — 2 spinodal line is 
absent. Contrarily, the v = 2 spinodal line based on extrapo- 
lating the dilute approximation to higher densities (see text) 
is shown by dash-dotted curves. 

The resulting phase diagram is shown in Fig. 2 for 
v = 2, 3, 4, 7, 10 spin states. The lower (T<T F ) and right 
(fci?|a|^,l) corner of the phase diagram is the spinodally 
unstable region where the system collapses. The region 
decreases for fewer spin states and is absent for v = 2. 
For comparison the spinodal lines for v = 2 and v = 3 
are shown when the dilute approximation of Eq. (|7|) is 
extrapolated into intermediate densities. Generally, the 
spinodally unstable regions based on the dilute approxi- 
mation are substantially overestimated. 

The regions of superfluidity given by T c = (j/tv)A and 
Eq. (|5|) are also shown in Fig. 2. As for the spinodally 
unstable region it is the lower right corner that is super- 
fluid. However, superfluidity extends to lower densities 
and therefore mechanical instability does not prevent the 
BCS-type pairing in the case of fermions. The opposite 
conclusion was reached for the pairing transition in Bosc- 
Einstein condensates pi]. As cooling becomes increas- 



ingly difficult at temperatures below the Fermi temper- 
ature we observe that the superfluid transition is readily 
obtained by increasing the density above fc^M^d- 

The phase diagram is quantitatively correct at low as 
well as high temperatures. Around the Fermi energy it 
gives a qualitative description only due to the approxi- 
mate thermal pressure employed. Furthermore, at inter- 
mediate densities the superfluid gaps become large ex- 
ceeding the Fermi energy for large spins, and the correc- 
tions to the ground state energies can therefore not be 
ignored. 



V. FINITE SYSTEMS 

The degenerate Fermi gases and Bose-Einstein conden- 
sates (BEC) produced so far contain n ~ 10 3 — 10 6 mag- 
netically trapped alkali atoms. Some of them interact via 
long scattering lengths such as the triplet 6 Li fermions 
with a = —2160 Bohr radii and singlet 85 i?&2 bosons with 
|a|^10 3 Bohr radii. Large scattering lengths a — > ±oo 
can be taylored by using Feshback resonances, i.e. hyper- 
fine states close to threshold further tuned by magnetic 
fields. Fermi gases differ from BEC's in several respects. 
Most importantly, whereas bosons sit at zero momentum 
states, fermions have considerable kinetic energy. There- 
fore, when interactions are small, a BEC has energy per 
particle Hlo and size a osc , where u> — (wi^) 1 ' 3 is the 
geometric average of the magnetic trap frequencies and 
a osc = h/^/muj is the oscillator length. A degenerate gas 
of N Fermionic atoms has a larger energy per particle 
~ N^-^huj and size L ~ N x / 6 a osc . In current experi- 
ments with degenerate Fermi gases and BEC's the den- 
sities are low so that the dilute potential applies and the 
energy per particle is approximately 
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(16) 



where the average density in the trap has been approx- 
imated by (p) ~ N/L 3 jL7). For a small number of 
trapped atoms with attractive scattering lengths the sys- 
tem is meta-stable but for a large number of trapped 
atoms, N^v{a osc / \v — l)a) 6 , the attractive potential 
overcomes the Fermi kinetic energy and the degenerate 
Fermi gas becomes unstable and implodes. However, 
around the same density /cp|a|£l and we enter the in- 
termediate density region, where the potential of (fill ) 
should be applied instead of the dilute potential. The 
gas is therefore mechanically stable for two hypcrspins 
only contrary to conclusions based on the dilute poten- 
tiallJU. 

Recent experiment on cold magnetically trapped 
Fermionic atoms observed egeneracy for 40 K atoms 
in the two hyperfine states tuf = 9/2, 7/2. Current ex- 
perimental oscillator lengths a osc ~ pm are less than 
one order of magnitude longer than the atomic scatter- 
ing length \a\ of 6 Li. It should be possible to reach in- 
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termediate densities, fcp|a|^,l, by trapping N^IO Li 
atoms rpTsjl . The atomic gases offer the unique opportu- 
nity to vary the densities as well as the relative amount 
of the hypcrfine states. Varying the composition is a con- 
venient wa y t o vary the gaps and attractive potential of 
Eqs. (f7|, p] , ^0| ) through v for given density and scattering 
length. In the limit where most atoms are in one of the 
states, the Fock and Hartree terms almost cancels and 
effectively v — » 1 + . 

More intricate systems of mixtures of fermions and 
bosons, e.g. 39 > 40 > 41 if isotopes can also be studied. If 
the interaction is attractive it will contract the atomic 
cloud towards higher densities. Irrespective of whether 
the bosons or fermions attract or repel the induced inter- 
actions, which are of second order in the fermion-boson 
coupling, enhance the gap ||. 

An artificial "gravitational" or "Coulomb" force can 
be exerted on the atoms by shining laser light on the 
trapped cloud from many directions |19(| . It would add 
an energy per particle of order ~ GNrr? /L to Eq. ( |l6| ) 
where G is proportional to the laser field intensity. Such 
an interaction has several interesting consequences. If G 
is attractive, it would contract the cloud towards higher 
densities which would increase gaps (see also |Q) and 
for sufficiently large G the intermediate density region 
is entered. Depending on the strengths and sign of the 
scattering amplitude and gravitational interactions the 
kinetic energy of the atoms will be balanced by the mag- 
netic trap and/or the scattering or gravitational interac- 
tions. The resulting phase diagram is much more com- 
plex. 

If such a strong attractive laser field is suddenly ap- 
plied to the gas, the Jeans instability sets in and the 
gas collapses until balanced again by the kinetic energy. 
Subsequently, the system will "bounce" analogous to the 
initial stages of a supernova explosion. If, however, in- 
termediate energies are reached and the number of spins 
exceed ^2.5, then the collapse will be further acceler- 
ated by the attraction between atoms. The correspond- 
ing critical particle number is 



N c ~ {Gm 3 a)- 3/2 



(17) 



at zero temperature. It differs from the standard Chan- 
drasekhar mass by a factor (m|a|) 3//2 because the insta- 
bility condition is kp\a\ ~ 1 whereas stars go unstable 
when the particles become relativistic kp ~ m. 



VI. SUMMARY 

The energy per particle and superfluid gaps have been 
calculated for an homogeneous system of fermions inter- 
acting via a long attractive s-wave scattering length. In 
the intermediate region of densities, where the interparti- 
cle spacing (~ 1/kp) is much longer than the range of the 
interaction but much shorter than the scattering length 
or |a|, the energy per particle and superfluid gaps are 



proportional to the Fermi energy. The energy per parti- 
cle increases linearly with the spin-isospin or hyperspin 
statistical factor such that, e.g., symmetric nuclear mat- 
ter is unstable in the intermediate density regions and 
undergoes spinodal decomposition whereas neutron mat- 
ter and Fermionic atomic gases with few hyperspin states 
are mechanically stable. 

A phase diagram of Fermi gases at low and intermedi- 
ate densities was constructed by including thermal pres- 
sures in the high temperature classical approximation. 
With the proper energy per particle at intermediate den- 
sities the spinodal region in the phase diagram was re- 
duced substantially as compared to conjectures based on 
extrapolations from the dilute limit. Generally, mechan- 
ical instability does not prevent a superfluid transition 
for a wide range of densities. This is contrary to Bose 
gases, where spinodal instabilities exclude pairing tran- 
sitions pl[ . 

The interaction energies of the many-body system were 
discussed for magnetically trapped cold degenerate gases 
of Fermi atoms. In such systems both superfluidity and 
the intermediate density region should be attainable. In 
these "novel" density regions the superfluid gaps can be 
large and the stability and sensitivity to the statistical 
factor v can be studied. Adding a gravitationally like 
force by shining laser light on the atomic cloud further 
increase densities whereby collapse and bounce analogous 
to the early stages of supernova explosions may be stud- 
ied. 
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